Neutrinos: the Key to UHE Cosmic Rays 
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Observations of ultrahigh energy cosmic rays (UHECR) do not uniquely determine both the 
injection spectrum and the evolution model for UHECR sources - primarily because interactions 
during propagation obscure the early Universe from direct observation. Detection of neutrinos 
produced in those same interactions, coupled with UHECR results, would provide a full description 
of UHECR source properties. 



PACS numbers: 98.70.Sa, 13.85Tp, 98.80.Es 
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Ultrahigh energy cosmic rays are one of the ten mir- 
acles of today's physics. Highest energy particles have 
long been suspected |lj to be extragalactic, because our 
Galaxy can not magnetically contain, and respectively 
accelerate, them. Current experimental data on the high- 
est energy cosmic rays show two severe problems: 

1. It is equally difficult to explain their production with 
either traditional astrophysical acceleration models [2J or 
with exotic top-down [3] particle physics models. Accel- 
eration models predict that UHECR are charged nuclei, 
whereas top-down models predict them to be 7-rays and 
neutrinos. 

2. Both nuclei and 7-rays have small energy loss dis- 
tances Ci oss = (l/E)dE/dx. Protons of 3xl0 20 eV lose 
their energy in propagation on 23 Mpc. Ci oss declines to 
less than 15 Mpc at higher energy. The 7-ray energy loss 
distance is less certain because of importance of the un- 
known isotropic radio background radiation. Reasonable 
estimates Q yield about 5 Mpc at energies between 10 19 
and 10 21 eV. Sources of UHECR then must be within 
several tens of Mpc, but none is identified. 

The solution of these problems is also affected by the 
current inconsistency in the results of the two major ex- 
perimental groups 0, E|- HiRes data seem to show a 
GZK cutoff, while AGASA's UHE cosmic ray spec- 
trum can be explained only with the addition of nearby 
sources or top-down scenarios. 

Observations of the UHECR energy spectrum do not 
uniquely determine the extragalactic cosmic ray source 
distribution or the source spectrum even with high statis- 
tics - there are too many different ways to fit the spec- 
trum. Figure ^ shows two extreme fits. The top panel 
illustrates a fit with a flat E~ 2 injection spectrum and 
(1 + z) m (m=3,4) cosmological evolution of the cosmic 
ray sources. Such fits 0, || can represent UHECR spec- 
trum adequately if galactic cosmic rays extend above 10 19 
eV. Correspondingly, the chemical composition of cosmic 
rays contains a fraction of heavy nuclei up to that energy. 
The second cosmic ray knee is where extragalactic 
cosmic rays prevail over the galactic ones. Fits with flat 
injection spectra require some cosmological evolution of 
the UHECR sources. Note, however, the small effect of 
cosmological evolution on UHECR above 10 18 eV - large 
redshifts do not contribute much. 

A different type of fit is illustrated in the bottom 
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FIG. 1: a) A fit of UHECR spectrum with flat extragalactic 
injection spectrum. Upper edge of shaded area is for m—4, 
lower - m=3. Galactic contribution is (dotted line) with 
7=2.6. b) Fitting with steep injection spectrum (7=1.7). 

panel of Fig. ^ Extragalactic cosmic rays with injection 
spectrum E -2 7 prevail down to 10 18 eV. These cosmic 
rays are expected to be protons and He nuclei. Gen- 
erally injection spectra steeper than E~ 2 - 5 can fit the 
observed spectra without strong redshift evolution. The 
second knee is the result of dE/dx due to pair produc- 
tion 0, 0| . The galactic cosmic ray spectrum extends 
to about 10 18 eV. The spectral shape of the extragalac- 
tic cosmic rays has to be flatter below 10 18 eV to avoid 
overproducing lower energy cosmic rays. One possible 
explanation is the limited horizon of lower energy extra- 
galactic cosmic rays because of scattering in extragalactic 
magnetic fields 
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Cosmogenic ( GZK) neutrino production. As shown in 
Fig. different UHE cosmic ray source models, with 
different luminosities at 10 19 eV, produce spectra consis- 
tent with observations. The differences at lower energy 
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function of red shift, we model the UHECR injection rate 
as 
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FIG. 2: Neutrino production in three simplified models for 
the source spectrum and evolution of UHECR. To enhance 
the contrast between different models, the source luminosity 
is extended out to 1 + z max = 10. 



may be disguised by contributions from the galaxy, or 
by propagation effects. The model degeneracy can be 
broken by considering the neutrino flux produced by the 
UHECR via their interactions on the cosmic micro-wave 
background |14Ul5l ]. the so-called GZK neutrino flux. 

To make the point, we consider a set of simple power 
law models in a matter dominated cosmology, with a ho- 
mogeneous distribution of UHECR sources. On cosmo- 
logical time scales, the UHECR interact quickly with the 
CMB, so we model the production of GZK neutrinos as 
instantaneous with particle injection. We use the neu- 
trino yields of Ref . [16j . The neutrino flux at Earth due 
to GZK production is 

d<£> f dT dv 

^E v ) = jdtde p -E v ^-(E v ,e p ,t) (1) 

where T is the injection rate of UHECR and y is the 
neutrino yield per proton injected with energy e p , and 
E v is the neutrino energy today. 

Equation ^ can be put into a more convenient form 
by defining q = 1 + z and integrating over redshift. For 
a matter dominated cosmology, H$dt = —q~ 3 / 2 d(\nq), 
where Hq is the present day Hubble parameter. As a 
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where 7 is the integral spectral index of the source, m 
describes the evolution of the co-moving source density, 
and A is chosen to normalize the emissivity of UHECR 
sources to their energy density in the present Universe. 
We use C-cr = 4.5 x 10 44 erg/Mpc 3 /yr as estimated by 
Waxman [TtJ • The neutrino yield can also be scaled with 
redshift. Defining the present day yield as ^-(E v ,E p ), 
the yield from a previous epoch is 



Eu-j^r{E u ,e p , z) 



E„^-(q 2 E v ,qe p ) (3) 



The factors of q are due to the redshift of neutrino en- 
ergy from production to the current epoch, and to the 
lowering of the reaction threshold due to the increased 
CMB temperature in the early Universe. Let Y$ be the 
integrated yield from injecting an E p < ' 1+1 ' > spectrum to- 
day: E U ^{E V ) = J 'dE p E p {1+l) E u ^{E„,E p ). Then 
the integral yield at any other redshift is given by 

J de p e-^E v -^-(E u ,e p ,t)=q^E u ^l( q 2 E u ) (4) 

With these definitions, the GZK production integral can 
be expressed as an integral over redshift, 



E V — {E V ) = — 

dE v Hq jq 



dE u 



(5) 

Written as an integral over In q, it is straightforward to 
see which epoch dominates the neutrino flux. If m + 7 = 
1.5 then all redshift intervals contribute with compara- 
ble importance. This is illustrated by the top panel of 
Fig. where m = 0.5,7 — 1- The thin lines represent 
the contribution to the neutrino flux from epochs spaced 
equally in lng, and of equal width in d(h\q). The peak 
contribution for each interval occurs at an energy E p k (q) , 
which scales with redshift as E p k(q) = E p k(Q)/q 2 . The 
sum of the thin lines, back to a red shift of l + z rnax = 10, 
gives the dotted curve - the predicted GZK neutrino flux 
for the model. The integrated flux is flat, with a width 
in E v of order (1 + z max ) 2 . 

Similarly, the middle panel illustrates a second model 
with equal contribution per epoch, although here there 
is no evolution (m = 0) and the spectrum is correspond- 
ingly steeper (7 = 1.5). The total cosmogenic neutrino 
flux is slightly lower because of the smaller number of 
interacting protons. The larger value of 7 is only evident 
through the different slope of the high energy part of the 
neutrino spectrum. In contrast, if m + 7 > 1.5 the neu- 
trino flux is dominated by past epochs, as illustrated in 
the bottom panel of Fig. [21 This corresponds to the sit- 
uation for a flat source model with significant evolution, 
to = 3,7 = 1. Note that this illustrative model is not 
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FIG. 3: Neutrino production for the UHECR models shown in 
Fig. The injection spectrum is cutoff at E c = 3 x 10 21 GeV 
and the source evolution breaks at z max = 3. Model results 
depend mildly on the cosmological model, = 0.0 or 0.7. 
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Discussion. The scaling analysis illustrates the combined 
importance of the UHECR spectrum and cosmological 
evolution in estimating the neutrino flux due to GZK 
production. For more realistic models, it is appropriate 
to consider other effects, such as a cosmological constant, 
a cutoff to the cosmic ray injection spectrum, and a cutoff 
or break in scale to the source evolution model. These 
effects alter the details of GZK production, but do not 
change the overall picture. Accordingly, Fig. [3] shows 
an order of magnitude difference in the GZK production 
for the two extreme UHECR models discussed in Fig. ^ 
Apart from the lower value of z max , the steep source 



model with no evolution (7 = 1.7, m = 0) is similar to 
the middle panel of Fig. [3 The model with flat spectrum 
and cosmological evolution (7 = 1, m = 3) predicts an 
order of magnitude larger flux. 

This brings us to the main point of this paper. The 
AUGER observatory, under construction, will measure 
the intensity and spectrum of UHECR. Above the GZK 
cutoff these are observations of the current Universe. Be- 
low the GZK cutoff it may be difficult to separate the 
old intergalactic cosmic rays from a population of young 
cosmic rays which originate within our own galaxy. Com- 
plementary to AUGER, experiments are being designed 
and constructed (ANITA, IceCube, Mediterranean km 3 , 
RICE), which may confirm the existence of GZK neutri- 
nos. A next generation of experiments (EUSO, OWL, 
SalSA, X-RICE) is being planned which would provide 
sufficient statistics (10-100 GZK events per yr) to comple- 
ment and expand the AUGER observations. Successful 
completion of one such experiment would be an impor- 
tant step toward understanding the sources of the highest 
energy particles in the Universe. 

It is often argued that neutrino astronomy has value in 
that neutrinos allow observation of the interior of objects, 
whereas photons only allow observation of the surface. 
Solar neutrinos are an example, confirming theoretical 
models of the thermonuclear furnace in the sun, by direct 
observation of the core. The current discussion is similar. 
The "surface" is the interaction distance for super-GZK 
UHECR. The "interior" is the early Universe, where the 
evolution of UHECR sources can be directly observed 
through the GZK neutrino flux. 
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